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Abstract. In this note, we provide a non trivial example of differential equation driven 
by a fractional Brownian motion with Hurst parameter 1/3 < iJ < 1/2, whose solution 
admits a smooth density with respect to Lebesgue's measure. The result is obtained 
through the use of an explicit representation of the solution when the vector fields of the 
equation are nilpotent, plus a Norris type lemma in the rough paths context. 



1. Introduction 

Let B = {B^, . . . , B"^) be a d-dimensional fractional Brownian motion with Hurst pa- 
rameter 1/3 < H < 1/2, defined on a complete probability space J-", P). Remind that 
this means that all the component B^ of B are independent centered Gaussian processes 
with covariance 

Rnit, s) := E [Bl 5^] = + t'^ - \t - 8^). (1) 

In particular, the paths of B are 7-Holder continuous for all 7 G {0,H). This paper is 
concerned with a class of M^-valued stochastic differential equations driven by 5, of the 
form 

d 

dyt = dBl, t G [0, T], yo = a, (2) 

i=l 

where T > is a fixed time horizon, a G M"* stands for a given initial condition and 
{Vi, . . . , Vd) is a family of smooth vector fields of M*". 

Stochastic differential systems driven by fractional Brownian motion have been the 
object of intensive studies during the past decade, both for their theoretical interest and 
for the wide range of application they open, covering for instance finance [T5| l32] or 
biophysics [20} [29] situations. The first aim in the theory has thus been to settle some 
reasonable tools allowing to solve equations of type ([2]). This has been achieved, when 
the Hurst parameter H of the underlying fBm is > 1/2, thanks to methods of fractional 
integration [271 123], or simply by means of Young type integration (see e.g [H]). When 
one moves to more irregular cases, namely H < 1/2, the standard method by now in order 
to solve equations like ([2]) relies on rough paths considerations, as explained for instance 
in da HUES] . 

A second natural step in the study of fractional differential systems consists in estab- 
lishing some properties about their probability law. Some substitute for the semigroup 
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property governing C{yt) in the Markovian case (namely when H = 1/2) have been given 
in [21 [21], in terms of asymptotic expansions in a neighborhood of t = 0. nSome consider- 
able efforts have also been made in order to analyze the density of C{yt) with respect to 
Lebesgue measure. To that respect, in the regular case H > 1/2 the situation is rather 
clear: the existence of a density is shown in [28] under some standard nondegeneracy 
conditions, the smoothness of the density is established in [19] under elliptic conditions 
on the coefficients, and this result is extended to the hypoelliptic case in [5]. In all, this 
set of results replicates what has been obtained for the usual Brownian motion, at the 
price of huge technical complications. 

In the irregular case H < 1/2, the picture is far from being so complete. Indeed, 
the existence part of the density results have been thoroughly studied under elliptic and 
Hormander conditions (see [6l [T2] for a complete review). However, when one wishes to 
establish the smoothness of the density, some strong moment assumptions on the inverse 
of the Malliavin derivative of yt are usually required. These moment estimates are still 
an important open question in the field, as well as the smoothness of density for random 
variables like yt. 

The current paper proposes to make a step in this direction, and we wish to prove that 
C{yt) can be decomposed as pt{z) dz for a smooth function pt in some special non trivial 
examples of equation ([2]). Namely, we will handle in the sequel the case of nilpotent vector 
fields Vi, . . . , V"(i (see Hypothesis 14. II for a precise description), and in this context we shall 
derive the following density result: 

Theorem 1.1. Suppose that the vector fields Vi, 1 = l,2,...,d are smooth with all 
derivatives bounded, and that they are n-nilpotent in the sense that their Lie brackets 
of order n vanish for some positive integer n. We also assume that Vi, . . . , satisfy 
Hormander 's hypoelliptic condition (their Lie brackets generate M™ at any point x G W^), 
and that all the Lie brackets of order greater or equal to 2 are constant. Then for allt > 
the probability law of the random variable yt, defined by ^ admits a smooth density with 
respect to Lebesgue measure. 

Notice that the hypoelliptic assumption is quite natural in our context. Indeed, it would 
certainly be too restrictive to consider a family of vector fields Vi, . . . , being nilpotent 
and elliptic at the same time. Moreover, some interesting examples of equations satisfying 
our standing assumptions will be given below. It should be stressed however that the basic 
aim of this article is to prove that smoothness of density results can be obtained for rough 
differential equations driven by a fractional Brownian motion in some specific situations, 
even if the general hypoelliptic case is still an important open problem. We refer to ^ for 
another case, based on skew-symmetric properties, where a similar theorem holds true. 

In order to prove Theorem II. H two main ingredients have to be highlighted: 
(i) Working under the nilpotent assumptions described above enables to use a Strichartz 
type representation for the solution to our equation, given in terms of a finite chaos 
expansion. This allows to derive some bounds for the moments of both yt and its Malliavin 
derivative, which is the main missing tool on the way to smoothness of density results for 
rough differential equations in the general case. 

(a) With the integrability of Malliavin derivative in hand, we shall follow the standard 
probabilistic way to prove smoothness of density under Hormander 's conditions, for which 
we refer to [T6 | l23 | 125] . To this purpose, the second main ingredient is a Norris type lemma. 
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which has to be extended (in the rough path context) to controlled processes. It should be 
mentioned at this point that a similar result has been proven recently (an independently) 
in dH. 

These two ingredients will be developed in the remainder of the article. 

Here is how our article is structured: Some preliminaries on rough differential equations 
and fractional Brownian motion are given in Section |5J Section [3] is devoted to the proof 
of our Norris type lemma for controlled processes in the sense of [T3]. Finally, Malliavin 
calculus tools and their application to density results for the random variable yt are 
presented at Section HI 

Notation: In the remainder of the article, c, Ci, C2 will stand for generic positive constants 
which may change from line to line. We also write a < 6 (resp. a x fe) when a < cb (resp. 
a = cb) for a universal constant c. 

2. Rough differential equations and fractional Brownian motion 

Generalized integrals will be needed in the sequel in order to define and solve equations 
of the form ([2]), and also to get an equivalent of Norris lemma in our context. Though 
all those elements might be obtained within the landmark of usual rough paths setting 
[T2| [22] we have chosen here to work with the algebraic integration framework, which 
(from our point of view) is more amenable to handy calculations. 

In this section, we recall thus the main concepts of algebraic integration. Namely, we 
state the definition of the spaces of increments, of the operator 6, and its inverse called 
A (or sewing map). We also recall some elementary but useful algebraic relations on the 
spaces of increments. The interested reader is sent to [H] for a complete account on the 
topic, or to IHl [13] for a more detailed summary. 

2.1. Increments. The extended integral we deal with is based on the notion of incre- 
ments, together with an elementary operator 6 acting on them. 

The notion of increment can be introduced in the following way: for two arbitrary real 
numbers £2 > ii > 0, a. vector space V, and an integer k > 1, we denote by Ck{[ii,i2]', V) 
the set of continuous functions g : [^1,^2]'^ V such that Qti - t^ = whenever ti = ti+i 
for some i G {0, . . . ,k — 1}. Such a function will be called a (/c — 1) -increment., and we 
will set C*([£i,£2]; = Ufe>iCfc([£i, £2]; V)- To simplify the notation, we will write CkiV), 
if there is no ambiguity about [£i,i'2]- 

The operator 5 is an operator acting on /c-increments, and is defined as follows on 
Ck{V): 

k+l 

S : C,{V) ^ C,+,{V), {5g)t,..,,^, = J2^-iy9t,...u-t,^,, (3) 

1=1 

where t^ means that this particular argument is omitted. Then a fundamental property 
of 6, which is easily verified, is that 66 = 0, where 66 is considered as an operator from 
Ck{V) to Ck+2{,V). We will denote ZCk{V) = Ck{V) n Ker5 and BCk{V) = Ck{V) n lm6. 

Some simple examples of actions of 6, which will be the ones we will really use through- 
out the article, are obtained by letting g G CiiV) and h G C2(y). Then, for any 
t,u,s G [ii,i2], we have 

{6g)st = Qt- Qs and {6h)sut = hst - hsu - Kt- (4) 
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Our future discussions will mainly rely on /c-increments with k = 2 or k = 3, for which 
we will use some analytical assumptions. Namely, we measure the size of these increments 
by Holder norms defined in the following way: for / G C2(y) let 

11/11, = sup and C^{V) = {/ G C^iV); < oo} . (5) 

Using this notation, we define in a natural way C^iV) = {/ G CiiV); \\Sf\\^ < oo}. In 
the sequel, we also handle norms including supremums, of the form 

||/IU,oo = 11/11, + II/II00, and Cr'°(\/) = {/GCi(\/); ||/||,,oo<oo}. (6) 
In the same way, for /i G C3(y) we set 



h\\^^p = sup 

s,u,te[ii,£2] 





hsut 




\U - S|T t 


— u\ 


\p 



(7) 



\\h\\^ = inf I ^ II /ii II p„;.-pj /i = ^ /li, < < j , 

where the last infimum is taken over all sequences {hi, z G N} C CslV) such that h = hi 
and over all choices of the numbers pi G (0,/i). Then ||-||, is easily seen to be a norm on 
Cs{V), and we define 

Cl^iV) := {heC,{V); \\h\\,<oo}. 

Eventually, let Cl'^(V) = U^^iCi^iV), and note that the same kind of norms can be 
considered on the spaces ZC^lV), leading to the definition of the spaces ZCt^{V) and 
ZCl~^{V). In order to avoid ambiguities, we sometimes denote in the following by J\f[ ■ ; Cj] 
the K-Holder norm on the space Cj, for j = 1, 2, 3. For ( G CjiV), we also set Af[C; CjiV)] = 

SUPsel£i;i2]'\\Cs\\v- 

The invertibility of 6 under Holder regularity conditions is an essential tool for the 
construction of our generalized integrals, and can be summarized as follows: 

Theorem 2.1 (The sewing map). Let fi > 1. For any h G ZCl^iV), there exists a unique 
Ah G Cl^iV) such that 6{Ah) = h. Furthermore, 

\m,<j^,m;C!;{v)]. (8) 

This gives rise to a continuous linear map A : ZC^iV) — )■ Ci^iV) such that 5 A = id^ct^iy)- 

Proof. The original proof of this result can be found in [13]. We refer to [HI [13] for two 
simplified versions. 

□ 

The sewing map creates a first link between the structures we just introduced and the 
problem of integration of irregular functions: 

Corollary 2.2 (Integration of small increments). For any 1-increment g G C2{V) such 
that 6g G Cg"*", set h = {id — A6)g. Then, there exists f G Ci{V) such that h = 6f and 

n 
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where the limit is over any partition Ust = {tQ = s, . . . ,tn = t} of [s, t] whose mesh tends 
to zero. The 1-increment 6f is the indefinite integral of the 1-increment g. 

We also need some product rules for the operator 6. For this recall the following 
convention: for g E C„([^i, £2]^'''^) and h E Cm([£i, £2]; M"^'^) let gh be the element of 
C„+„_i([£i,£2];K''^) defined by 

{gh)tl,...,tm + n-l = gtl,...,tnht„,...,t^ + „-.l (9) 

for ti, . . . ,tm+n-i £ [■^i,-^2]- With this notation, the following elementary rule holds true: 

Proposition 2.3. Let g E C2{[ii, £2];'^^''^) and h E Ci{[ii, £2]]'^'^) ■ Then gh is an element 
o/C2([£i,^2];K0 and5{gh) =5gh-g5h. 

2.2. Random differential equations. One of the main appeals of the algebraic inte- 
gration theory is that differential equations driven by a 7-Holder signal x can be defined 
and solved rather quickly in this setting. In the case of an Holder exponent 7 > 1/3, the 
required structures are just the notion of controlled processes and the Levy area based on 

X. 

Indeed, recall that we wish to consider an equation of the form 

d 

dyt = J2^Mdxi, tE[0,T], yo = a, (10) 

i=l 

where a is a given initial condition in M"^, x is an element of C7([0, T]; M"^), and {Vi, . . . , Vd) 
is a family of smooth vector fields of M"^. Then it is natural that the increments of 
a candidate for a solution to ( llOp should be controlled by the increments of x in the 
following way: 

Definition 2.4. Let z be a path in (M*") with 1/3 < k < 7, and set 5x := x-*-. We say 
that z is a weakly controlled path based on x if Zq = a with a E M"^, and 6z E (M"*) has 
a decomposition 6z = (x.^ + r, that is, for any s,t E [0, T], 

6zst = Ci^f +rst, (11) 

where we have used the summation over repeated indices convention, and with C"*^, . . . , G 
(M"'), as well as r E C|^(M™). 

The space of weakly controlled paths will be denoted by Q^^(M'"), and a process z E 
Q^^(M"*) can be considered in fact as a couple {z,Q. The space Q^ai^"^) endowed 
with a natural semi-norm given by 

d 

Af[z; Ql^n] = Ar[z;C^iRn] + J2Af[C;Ct'\Rn] + Ar[r;Ci''{Rni (12) 

i=i 

where the quantities JVlg; Cj] have been defined in Section [2771 For the Levy area associ- 
ated to X we assume the following structure: 

Hypothesis 2.5. The path x : [0, T] — M'^ is ^-Holder continuous with | < 7 < 1 and 
admits a so-called Levy area, that is, a process x^ G C2^{M.'^''^), which satisfies 5x^ = 
x-*- ® x-*-, namely 

for any s,u,t E [0, T] and i,j E {1, . . . ,d}. 
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To illustrate the idea behind the construction of the generalized integral assume that 
the paths x and z are smooth and also for simplicity that d = m = 1. Then the Riemann- 
Stieltjes integral of z with respect to x is well defined and we have 

/"* 1 

I Z^dXu — Zg{xt Xg^ -\- I (^Z^ Zg^dXu — ^s'^st ~^ / i^^^^sudXu 
J s J s J s 

for £i < s < t < £2- If ^ admits the decomposition (ITT]) we obtain 

/ {5z)sudxu= I {Cs^lu + Psu) dxu = Cs x.l^dxu+ / psudXu- (13) 

Js J s J s J s 



Moreover, if we set 



^st 



< s<t< 



then it is quickly verified that x is the Levy area associated to x. Hence we can write 



/"* 12/"* 

/ Zqj,dXu ^s^sz ~^ Cs ^st ~^ / Psu dXu- 
J s J s 



Now recast this equation as 



/"* /"* 1 

/ psu d^u I '^u ZgX.^^ 
J s J s 



(14) 



and apply the increment operator 5 to both sides of this equation. For smooth paths z 
and X we have 



5 



z dx 



0, 



5{z^^] 



-52; x"*", 



by Proposition 12.31 (recall also our convention (|9]) on products of increments). Hence 
applying these relations to the right hand side of (fT4|) . using the decomposition ( ITT]) , the 
properties of the Levy area and again Proposition 12. 3[ we obtain 



p dx 



c2 

sut 



sut 



5zsu^lt + 5Csuy^lt-Cs5^ 



ut 



— Psu^ut ~^ ^Csu ^uf 

Summarizing, we have derived the representation 



pdx 



Psu^ut "I" ^Csu ^uf 



sut 



As we are dealing with smooth paths we have 5 (J pdx^ lies into the space ZC^'^ and thus 
belongs to the domain of A due to Proposition 12.11 (Recall that 66 = 0.) Hence, it follows 

/ psudxu = Ast {p^^ + 6Cx^) , 

J s 

and inserting this identity into (fT3|) we end up with 

Zudxu = Zs y.]t + Cs ^It + J^st {p + 6C x^) . 
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Since in addition 

px^ + 5Cx^ = -(5(2x^ + Cx2), 

we can also write this as 

Thus we have expressed the Riemann-Stieltjes integral of z with respect to x in terms of 
the sewing map A, the couple (x-*^, x^) and of increments of z. This can now be generalized 
to the non-smooth case. Note that Corollary 12 . 21 justifies the use of the notion of integral. 

Proposition 2.6. For fixed ^ < k < let x be a path satisfying Hypothesis \2.5\ on an 
arbitrary interval [0,T]. Furthermore, let z G ^{[ii, £2]',^.'^) such that the increments 
of z are given by / fli]) . Define z by % = a with d G M and 

Szst = [(id-A5)(z*xi'* + e^'^')],, (15) 

for ii < s < t < £2- Then J'{z* dx) := z is a well-defined element of q,([^i, £2]] ^) o-nd 
coincides with the usual Riemann integral, whenever z and x are smooth functions. 

Moreover, the Holder norm of J{z* dx) can be estimated in terms of the Holder norm 
of the integrator z. (For this and also for a proof of the above Proposition, see e.g. |14|.) 
This allows to use a fixed point argument to obtain the existence of a unique solution for 
rough differential equations. 

Theorem 2.7. For fixed ^ < k < 'j, let x be a path satisfying Hypothesis \2.5\ on an 
arbitrary interval [0,T]. Consider a given initial condition a in and (Vi,...,Vrf) 
a family of vector fields of , bounded with bounded derivatives. Let ||/||^,oo = 
ll/lloo + ll^/llft ^6 ^he usual Holder norm of a path f G Ci([0,T]; M}). Then we have: 

(1) Equation [W\) admits a unique solution y in ^([0, T]; M™) for any T > 0, and 
there exists a polynomial : — > M"*" such that 

^f[y;Q:AM;Rn] < Pt{\\^%A\^'\M (16) 

holds. 

(2) Let P : M™ X CJ{[0,T];R'^) x C2^([0, T]; M™'"^) CJ{[0,T];R"') be the mapping 
defined by 

F (a,x^x^) = y, 

where y is the unique solution of equation / flOj) . This mapping is locally Lipschitz 
continuous in the following sense: Let x be another driving rough path with corre- 
sponding Levy area x^ and a be another initial condition. Moreover denote by y 
the unique solution of the corresponding differential equation. Then, there exists 
an increasing function Kt '■ R^ — R'^ such that 

\\y - y\U,oo,T < /iT(||x^||^, llx^ll^, ||X^||27, 11X^1127) (17) 

X (|a — a| + llx"*" — x"*"!!^ + ||x^ — x^||27) 

holds. 

The theorem above is borrowed from [T2| [HI [22], and we send the reader to these 
references for more details on the topic. 
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2.3. Fractional Brownian motion. We shall recall here how the abstract Theorem 12.71 
applies to fractional Brownian motion. We will also give some basic notions on stochastic 
analysis with respect to fBm, mainly borrowed from [26], which will turn out to be useful 
in the sequel. 

As already mentioned in the introduction, on a finite interval [0, T] and for some fixed 
H e (1/3,1/2), we consider (fi,^-', P) the canonical probability space associated with 
fractional Brownian motion with Hurst parameter H. That is, Vt = Co([0, T]; M'^) is the 
Banach space of continuous functions vanishing at equipped with the supremum norm, 
J-" is the Borel sigma-algebra and P is the unique probability measure on VL such that 
the canonical process B = {Bt, t G [0,T]} is a rf-dimensional fractional Brownian motion 
with Hurst parameter H. Specifically, B has d independent coordinates, each one being 
a centered Gaussian process with covariance given by ([1]). 

2.3.1. Functional spaces. Let £ be the set of the space of d- dimensional elementary func- 
tions on [0, T]: 

rij — l 
i=0 

for j = l,...,rf}. (18) 
We call 7i the completion of S with respect to the semi-inner product 

d 

{fi 9)n = '^(fi^ 9t)Hi^ where (l[o,t], l[o,s])w, := ^(s, i), s,tG[0,T]. 
1=1 

Then, one constructs an isometry K'}^ : H /.^([0, 1]; M*^) such that 
where the kernel Kh is given by 

+ Q - iJ^ ^ v^'-^v - uf--^ dvj l{o<«<o, (19) 

with a strictly positive constant Ch, whose exact value is irrelevant for our purpose. 
Notice that this kernel verifies i?/^(t, s) = J^^^ KH{t,r)KH(s,r) dr. Moreover, observe 
that can be represented in the following form: for if = [ipi, . . . ,ipd) G Ti, we have 
K},^={K*j,^\...,K},^''),with 

K*hV={K*h^\---,K*h^'), where [K*h^%, = dnt'/'-'' [d]!^"'' {u~^'/'~''^^') 

for a strictly positive constant dn- In particular, each Tii is a fractional integral space of 
the formX^/'-^(L2([0,T])) and Cy'"^([0,T]) C Ki. 
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2.3.2. Malliavin derivatives. Let us start by defining the Wiener integral with respect to 
B: for any element f in £ whose expression is given as in ( 1T8|) . we define the Wiener 
integral of / with respect to B as 

d nj-l 
3=\ j=0 

We also denote this integral as J^^ f{t)dBt, since it coincides with a pathwise integral with 
respect to B. 

For 6* : M — 7- M, and j E {1, . . . ,d}, denote by 9^^^ the function with values in M°' having 
all the coordinates equal to zero except the j-th coordinate that equals to 6. It is readily 
seen that 



E 



Bil^ABhf^ 



[0,s) ^ V [0.*) 



This definition can be extended by linearity and closure to elements of "H, and we obtain 
the relation 

E[B{f)B{g)] = {f,gU, 
valid for any couple (/, g) G 'H?. In particular, B{-) defines an isometric map from "H into 
a subspace of L?'{VL). It should be pointed out that P) defines an abstract Wiener 

space, on which chaos decompositions can be settled. We do not develop this aspect of 
the theory for sake of conciseness, but we will use later the fact that all norms are 
equivalent on finite chaos. 

We can now proceed to the definition of Malliavin derivatives, for which we need an 
additional notation: 

Notation 2.8. For n,p > I, a function f E C^(]R";]R) and any tuple (ii,...ip) E 
{1, . . . , dy, we set di,...^J for 

With this notation in hand, let S be the family of smooth functionals F of the form 

F = /(i?(/ii),...,i?(/i„)), (20) 

where E Ti, n > 1, and / is a smooth function with polynomial growth, 

together with all its derivatives. Then, the Malliavin derivative of such a functional F is 
the "H-valued random variable defined by 

n 

DF = Y,d^f{B{hr), . . . ,B{K))K 
1=1 

For all p > 1, it is known that the operator D is closable from L^^Q) into U'iVt^'H) (see 
e.g. Chapter 1]). We will still denote by D the closure of this operator, whose domain 
is usually denoted by D^'^ and is defined as the completion of S with respect to the norm 

\\Fh,v ■= mF\n + E{\\DFrn))'^. 

It should also be noticed that partial Malliavin derivatives with respect to each component 
B^ of B will be invoked: they are defined, for a functional F of the form (120|) and 
j = 1, . . . as 

n 

D^F = Y,^^f{B{h,),...,B{K))hf, 

i=l 
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and then extended by closure arguments again. We refer to Chapter 1] again for the 
definition of higher derivatives and Sobolev spaces D'^'^ for k > 1. 

2.3.3. Levy area of fBm. There are many ways to define the Levy area associated to 
fBm, and the reader is referred to [T2l Chapter 15] for a complete review of these. The 
recent paper jTT] is however of special interest for us, since it enables a direct definition 
of by Wiener chaos techniques. It can be summarized in the following way: 

Proposition 2.9. Let 1/3 < H < 1/2 be a fixed Hurst parameter. Then the fBm B 
belongs almost surely to any space Cj for ^ < H , and it gives raise to an increment 
B^ G Cg'^ which satisfies Hypothesis \2.5[ Furthermore, for any < s < t < T , B^^ is an 
element of the second chaos associated to B, and 

E[|B?,n <c,(t-.r^ p>i. 

Moreover, the iterated integrals of B can be obtained as limits of Riemann type inte- 
grals. Indeed, for /c > 1 and < s < t < T, consider the simplex 

Sk{[s, t\) = {{ui, . . . ,Uk); s < ui < ■ ■ ■ < Uk < t} . (21) 

For a given partition 11 of [0, T], we also denote by B^ the linearization of B based on 11. 
Combining the results of [TT| [T^ . the following proposition holds true: 

Proposition 2.10. Let k >1, and for a sequence of partitions (n„)„>i, set B"- := B^" . 
For < s < t < T and {ii, . . . , ifc) G {1, . . . , d}^, we consider then 

gk n,n,...,i, ^ f . . . 

St I Ui Uf^ ' 

JSk{[s,t]) 

understood in the Riemann sense. Then there exists a sequence of partitions (n„ )„>i such 
that B'*^'"'*^'"''*'' converges almost surely and in , as an element ofd^"' for any < H, to 
an element called B'''*^'"''** . When k = 1, we obtain the increment 6B of our fBm. When 
k = 2, the limit corresponds to the increment B^ of Proposition IKVi 

As a corollary of the previous considerations, we have the 

Proposition 2.11. Assume 1/3 < H < 1/2. Then Theorem \2. 71 applies almost surely to 
the fBm paths, enhanced with the Levy area B^ . We are thus able to solve equation 

d 

dyt = Y,V,{yt)dBi, tG[0,T], yo = a, (22) 

i=l 

under the conditions of Theorem \2. 71 

3. A NORRIS TYPE LEMMA 

Norris' lemma [25] is one of the basic ingredients in order to obtain smoothness of 
densities for solutions to stochastic differential equations under hypoelliptic conditions, 
and was already extended to fBm with Hurst parameter i7 > 1/2 in [3]. We shall extend 
in the current section this lemma to the rough paths context. A preliminary step along 
this direction consists of proving the following elementary lemma: 



NILPOTENT RDES 



11 



Lemma 3.1. Let < a < p < 1, and consider b G C([0, T]). Then for any < r] < 1, we 

have 

\\b\\a,oo<C^,p[v\\bl,oo + V-'^^'~"^\\b\\L,] , (23) 
where we recall that ||6||q,oo has been defined by 

Proof. Recall that ||fe||Q, has been defined by ([5]). Thus, for < a < p < 1, we have 



sup 

s,t 



\Sb 



St 



p 



<2'-7\ 



-1 cx_ 

|oo"||6||;, 



7||6||p + 7 



\t - s\P 

Thus, for an arbitrary constant rj > 0, we have 

II^IU < Ca,p {v\\b\\p + \\b\\oo) < C^,p {v\\b\\p + V 

thanks to Young's inequality. Therefore for an arbitrary constant < ?7 < 1 

\\b\\a<C^,p{v\\b\\p + V-^\\b\\oo) . 
Invoke now the interpolation inequality [3], formula (3-17)] in order to get 

\\b\\a < C^,p{v\\b\\p + V~^- 

= C^AvWbWp + V-'^^'-'-^lbWL,] , 

where we have chosen 7 = ?7~. Invoking again [31 formula (3.17)] in order to go from 
II ■ II Q II ■ lU.oo norms and this finishes our proof. 

□ 

We can now turn to the announced Norris type lemma, whose proof is an adaptation 
of |3] to the case of controlled processes. 

Proposition 3.2. Assume B is a fractional Brownian motion with H > 1/3. Let z be a 

controlled path in Q!^{M.'^), with decomposition 

6zl, = Q^Blf+rl,. (24) 



We assume that 1/3 < a < 'y < H and that the quantity 'E[J\fP[z; Q:^(M'^)]] is finite for 
all p > 1. Set 6yst = J{z*dB) according to Proposition \2.b\ Then there exists g > such 
that, for every p > we can find a strictly positive constant Cp such that 

P (|||/|l7,oo ^ ^1 and ||^||q.,oo 

> 5") < CpeP. 

Proof. In order to avoid cumbersome indices, we shall prove our result in the case of 
1-dimensional processes. Generalization to the multidimensional setting is a matter of 
trivial considerations. We also work on the interval [0,1] instead of [0,T] for the sake 
of notational simplicity. As a last preliminary observation, note that if z admits the 
decomposition flM|) . then according to f fTSj) we have 

6y = zB^ + CB^ + y\ where = A (r + 5C B^) . 

Similarly to what is done in j3], we consider two time scales 5 ^ A ^ 1. We assume 
moreover that A/6 = r with r G N. We use a partition n < 1/5} of [0,1] with 
tn = n6, so that t^r = NA. Some increments below will then be frozen on the time scale 
A, in order to take advantage of some averaging properties of the process B. 
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Stepl: Coarse graining on increments. Consider then n such that (A^ — l)r <n< Nr — 1, 
so that (A^ - 1)A <tn< NA - 1. According to ([15]), we have 

^yt„t„+i = Bj^j^^^ + (t„ ^t„t„+i + ytntn+i 

= zma Bf , — 6zf MA B?" , + Ct B? , + Vi i , , 
where is an increment in €2^^ . Thus 

^^TVA ^l„u+i = ^yt,^t„+i + 52;t„,7VA ^tt„+i - Ct„ Bf^i^^^ - yl^tn+i ■ (25) 
Set now Xn = J2nL'{N-i)r l|B?„t„^J|^ and Yn = xH'^. Then 

\zna\^\Xn\= Y1 \^N^'^t„tn+i\ ■ 
n={N~l)r 

Furthermore, invoking relation fl25|) . we get 

Raising this inequahty to power 4 and summing over (A^ — l)r < n < Nr — 1, we obtain 

I4a Xn\ < S^-'-'A {\\y\\, + ||2||,||5||,A^ + llClloollB^lb,^^ + \\y%,6'^)' , 
and therefore 

kiVAllAT < S^-'/'A'/' {\\y\\, + ||2||,||5||,A^ + IIClloollB^lh,^^ + \\y%,S'^) • 
Sum now over A^ (recall that 1 < A^ < 1/A) in order to get 

l/A 

J2 I^^All^iv < (5^-i/'^A-3/4 (\\y\\^ + ||z||^||i?||^A^ + llClloollB^lls-,^^ + \\y%,S'^) . (26) 

N=l 

Step 2: Behavior of a order variation. Throughout the proof, we shall use the notations 
<, X given in the introduction. For K > 1^ set Xk = ^^=1 IB^^^^^J^. We shall prove 
that 

^[Xk] X K5^", and Var(X^) < KS"^^ . (27) 

Indeed, a simple scaling argument shows that Xk ^— Xk, with Xk = X]^=i l-^n n+il^ 
Introduce now the k-th Hermite polynomial (see [22] for a definition and properties 
of these objects), and notice in particular that H2{x) = — 1 and H^lx) = — 6x^ + 3. 
This enables to decompose Xk as 

K 

Xk = Y1 [H4{Bl,^+i) + 6H2{Bl^^,)] + 3K. (28) 

n=l 

Recall now that for a centered Gaussian vector (f/, V) in such that E[t/^] = E[y^] = 
1 we have 

E[Hk{U)] = 0, and E [Hk{U) Hi{V)] = k\ (ElUV])" ^k=iy (29) 
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= E[B^^,„^+, = - [\n2 - m + Ip^ + \n2 - ni - Ij^^ - 2|n2 - ml^^] . 



Plugging this identity in 0281) . this immediately yields E[X/^] = 3K, which is our first 
assertion in (127|) . In addition, the second part of ( 12^ entails 

K 

Var(X^) = 2 (12<„^ + <J :=2^^, 

ni,n2=l 

where we have set 

_ 1 

Summarizing, we have obtained that 

Var(XK) = S^" VsLTiXK) = 2 S^^ Sk- (30) 
We will now prove that Sk ^ K. Indeed, write first 

Sk= (12<,„, 2 (12 <,„, + <, J :=^], + 2^|. 

l<ni<ii: l<ni<n2<X 

Then, since am.ni = 1, it is readily checked that S]^ < K. Moreover, the term S'j^ can be 
decomposed into 

Sl= Y (l2<„,+i + <,„,+i)+ Y (12<,.2+<J :=^l'+^f- 

l<ni<A'— 1 l<ni,n2<K,n2—ni>2 

Notice now that ani,ni+i = ~[1 ~ 2^"^^^^^^"^^], which easily yields S]^ < CrK. 
As far as 5*^ is concerned, write for n2 — rii > 2 



pi pr 

H{2H-1) dr 

Jo J-1 



(n2 — ni) + u\^^ ^ du^ 



l4:H-4 



which immediately yields ani,n2 ^ 1(^2 — ''^i) — M"^^ ^- Thus 

Sk < E E \in2-n,)-ir-' + 3\in2-n,)-l\ 

ni=l n2=ni+2 

A'-2 AT-jii + l oo 
ni=l m=l m=l 

since m'^^"^ is finite whenever H < 3/A. 

Gathering our bounds on Sj^, S]^ and S]^, we obtain < K, and plugging this bound 
into (130|) . we end up with Var(X/^) < K 6^^ , which is our claim. This finishes the proof 
of ([27]). 

Step 3: Concentration inequalities for Y^- Let us recall that Xn is in the 4*^^ chaos of the 
fBm B. Hence, a result by Borell [5] entails 

|Xjv — E[Xjvj| \ -c9«i/2 - „ 

> n < cie ^ , M > 0, 



[Var(X;v)]'/' 

for two universal constant Ci,C2 > 0. With fl27|) in hand, this yields 

P (|Xjv - 3A54^-i| > Ai/2^'^/f-i/2^^) < cie-^^"''', u>0. (31) 
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We now wish to produce a concentration inequality for = . Since Xj^ is a small 
random quantity of order A^^"^"^, let us use the inequality 

|5i/4_^i/4| <_t-3/4|^_^|^ where ^ e {aAb, aWb). 
Apply this with ^ = lA^^"^"^ in order to get 

with 

A, = p(^|X^-3A5^^-i|>cA^/V^-i/V^iv>^A<5^^-i^ 

A2 = p(^Xj,<^-A6'^-'y 

Furthermore, a straightforward application of f l3T]) gives 

< Cle-^^"''^ and A, < c^e-^-^"-/'^'''. 

Plugging these inequalities into (132|) . we end up with the following concentration inequality 
for Yn. 

P {\Yn - 3V4aV45H-i/4| > c A-^/^5^+i/^m) < cie-'^^"''' + c.e-'''^^/'^'''. (33) 

We shall thus retain the fact that Y^ i a random quantity of order 3^/^A^^'^S^~^^^, with 
fluctuations of order A^^^'^S^'^^^^: 

\Yj^ - 3V4Al/45/^-l/4| ^ ^H+l/4^-l/4_ (34) 

Step 4- Use of the interpolation inequality. Start again from equation ( 126|) . One would 
like to have an approximation of the norm of z appearing on the left hand side of this 
inequality, that is one would like to replace Y/v by A. To this purpose, replace first Y^ by 
its approximation A^/^(5^~^/^ from the last step. This yields an inequality of the form: 

l/A 1/A 

J2 \zj,^\ < \\z\\^ J2 \Yn - 3V4Al/45^^-l/4| 
N=l N=l 

+ 5^-'/'A-3/'(||y||,+ lkll,||S||,A^ + ||CI|oo||B2||2,5^ + lb«||3,5'^). 
Rescale this inequality in order to get A multiplying on the left hand side, which gives: 

l/A 

A ^ \zna\ 

N=l 

< Iklloo RsA + '^"^''"^^ {hh + IWIBIU^"' + llClloollB^lh^^^ + \\y%,6^"') , (35) 

where 

l/A 

RsA ■■= (5-(^-i/^)A3/'i Yl - 3^/'Ai/^(5^-i/% (36) 

Af=l 
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Furthermore, it is well known that |A^]^^^ \zna\ ~ W^Wl'^I ^ II-2||7A'^7 so that we can 
recast (ESI) into 



\\z\\li < ll-zll^A^ + ||2||oo-R<5,A 

+ ^"^""'"^ {Wvh + \\44BIA' + ||Cl|oo||B^||2,5^ + \\y%,6''^) . (37) 

Recall that we take < a < 7 < if and set uh '■= 1/(7 — a). According to ( 125]) we 
have 

||^||a,oo ^ ^||2^||7,oo ~l~ ^ ll^llii) 

for any (small enough) constant 77, and plugging ( 137|) into this last relation we obtain 

||-2'||a,oo ^ ^7 A''' II 2^11 ^^Qo -|- Tj ||2^||oo -R<5,A ~l~ '7ll'^ll7,oo 
+ r'^"(lbll7,oor(^-^) + ||^||,,oo||i?||7r(^"^)A^+ ||C||oo||B^|h,5^^"^+ ||y»||375'^^''). 

Defining e as e = ii — 7, we get 

II^IU.oo ^ ''7 '^"A'^'ll^ll^^oo + ^7 ll^lloo -R<5,A + ^11^117,00 

+ r''"(lbll7,oor^>lkll7,oo||5||^r^A^+||C||oo||B2||2^5^-^>||y«||^ (38) 

Step 4- Tuning the parameters. Recall that we have chosen ^ 5 ^ A ^ 1. We express 
this fact in terms of powers of e, by taking 5 = and A = e^, with /i > A > 0. We 
also choose 77 of the form 77 = e^l^" . We shall now see how to choose A, /i, r conveniently: 
write (l38l) as 



Irll < r"^+^'^ll2ll +r"^ll2ll /?A A + r^/^^llzll + ll?/ll ^-^-^^ 

I II a, 00 rN^^ ||'''|l7,OOT^C ||<0||oo-' ''0,A TC ||'^||7,ooT||y II 7,ooC- 

\z\ 



„,,o||5|l7^-^-''^"''" + ||Cl|oo||B^||2,e-^+^(^-^") + |b«||37^-^"'^^'"-'^. (39) 

In order to be able to bound z when y is assumed to satisfy |||/||7,oo < ^1 the coefficients 
in the right hand side of (139|) should fulfill the following conditions: 

• The coefficient in front of Hyjl-y^oo should be smaller than . 

• The other coefficients should be ^ 1. 

Looking at the exponents in (1391) . assuming that e is arbitrarily small and letting for the 
moment Rs,a apart, this imposes the following relations: 

A7 > r, and < r < 1. (40) 

Let us go back now to the evaluation of Rs,a, given by expression (l36i) . with the order 
of magnitude of \Yn- — A^^^5^~^^^\ given by (IMl) . Therefore 

Expressing this in terms of powers of e, we end up with 

V'""" Rs,A ~ e'^, with K = ^^-^ - r. 
If we wish this remainder term to be small, this adds the condition 

/i - A > 2r, (41) 
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which can be fulfilled easily. From now on, we shall assume that both (HUj) and f HT]) are 
satisfied. 

Step 5: Conclusion. Recall that we wish to study the probability P(|||/||^^oo < and 
||-2||a,oo > This quantity is obviously bounded by Bi + B2, where 

B, = P {\\y\\^,^ < e, ||2|U,oo > \R5a\ < c A-i/^5^+i/^(A/<5)^") 

where e is an arbitrary small positive constant. Furthermore, inequalities (155]) and (15^ 
yield, for any p > 1, 

where we have used the fact that bj t^ = e^~^ . 

We can now bound B\. notice that according to fl39|) . if we are working on 

then there exists a p > such that 

lklU,oo < {\^M\z- Q^(M-)] + \\B% + llB^liy . 

Moreover, recall that M\z] Q:^(]R'")] is assumed to be an U random variable for all r > 1, 
while II II 7 and ||B^||2^ are also elements of U ^ since they can be bounded by a finite 
chaos random variable. Thus Tchebychev inequality can be applied here, which entails 

fii < (1 + E {M''\z; (M'")] + \\Bf^ + llB^ll^y ) 

for an arbitrary / > 1. It is now sufficient to choose q < p and / large enough so that 
l{p — q) = p to conclude the proof, by putting together our bounds on Bi and i?2. 

□ 

4. Malliavin calculus for solutions to fractional SDEs 

This section is the core of our paper, where we derive smoothness of density for the 
solution to f l22]) . We first recall some classical notions on representations of solutions to 
SDEs, and then move to Malliavin calculus considerations. 

4.1. Representation of solutions to SDEs. The first representations results for solu- 
tions to SDEs in terms of the driving vector fields can be traced back to the seminal work 
of Chen [7j. They have then been deeply analyzed in [T8l [30] . and also lye at the basis 
of the rough path theory [22]. We have chosen here to present these formulas according 
to [1], which is a recent and didactically useful account on the topic. 

Recall that we are considering a d-dimensional fBm B with 1/3 < H < 1/2. According 
to Section 12. 3[ this allows to construct some increments B'' out of B which can be seen 
as limits of Riemann iterated integrals over the simplex Sk{[s,t]), as recalled at Proposi- 
tion [5]Tni Furthermore, one can solve equation ( 12^ under the conditions of Theorem 12.71 

Let us introduce some additional notation: let V be the space of smooth bounded 
vector fields over M™. If G V, the vector field exp{V) G V is defined by the relation 
[exp(V)](^) = where t > 0} is solution to the ordinary differential equation 

dtMO = v{^m, ^o(o = e. (42) 
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The aim of Chen-Strichartz formula is to express the solution yt to equation f l22]) . for an 
arbitrary t G [0,T], as yt = [exp(Z()](a) for a certain Zt G V. 

To this purpose, let us give some more classical notations on vector fields: if G V, 
then the vector field [V, W] G V (called Lie bracket of V and W) is defined by 

Notice that this notion is usually introduced though the interpretation of V as a set of 
first order differential operators. A Lie bracket of order k can also be defined inductively 

for A; > 2 by setting 

[Ui---Uk] = [[U,---Uk^,],Uk], 

for Ui . . .Uk G V. With this notation in hand, our main assumption on the vector fields 
Vi, . . . ,Vd governing equation f l22|) is the following: 



Hypothesis 4.1. The vector fields Vi, . . . , are n-nilpotent for some given positive in- 
teger n. Namely, for any (ii, . . . , i„) G {1, . . . , c?}", we have [V^^ ■ ■ ■ l^,J = 0. 

We are now ready to state our formulation of Strichartz' identity, for which we need 
two last notations: for > 1, we call &k the set of permutations of {1, . . . , k}. Moreover, 
for cr G (5^, write e{a) for the quantity Card({j G {1, . . . , — 1}; cr(j) > cr(j + 1)}). Then 
the following formula is proven e.g. in [H [181 130] : 

Proposition 4.2. Under the hypothesis of Theorem \2. % let y be the solution to equa- 
tion (f£g|). Assume Hypothesis 4jJ_ holds true, and consider t G [0,T]. Then yt = 
[exp{Zt)]{a) , where Zt can be expressed as follows: 

k=l n,...,jfe = l ""GSfc \e(cr)/ 

where we have set r = and Vj^ j^. = [V^j ■ ■ ■ Vi^] in the formula above. 

As a warmup to the computations below, we prove now that one can extend our in- 
equality (fT6l) thanks to Strichartz representation, covering the case of unbounded vector 
fields with bounded derivatives: 



Proposition 4.3. Suppose Hypothesis \4.1\ holds true, and that the smooth vector fields 
Vi, i = 1,2, ■ ■ ■ ,d have bounded derivative. Assume moreover that all the Lie brackets of 
order greater or equal to 2 are bounded vector fields. Then the solution y of equation [2^] 
admits moments of any order. Namely, for any m > 1 and any T G (0, oo), 



E 



sup \ytr 

0<t<T 



< OO . (43) 



Proof. One can restate Proposition 14.21 as follows: for any t < T, the random variable yt 
can be expressed as yt = (pi, where 0* := 0^ : M"' — > M"' satisfies (for t fixed) 



s 

n—l d 



ds<Ps = Y. E v^r'-''^v,,„„,,(0,), o<s<i 
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Let us separate the first order integrals in this equation, in order to get 

d n—1 d 

i=l k=2 ii,...,ifc=l 

Since Vi, i = 1,2, ■ ■ ■ ,d, have bounded derivatives and since all the Lie brackets of order 
greater or equal to 2 are bounded we see that 

n-l 

\^ U/,n,-.,ifc|lA/. , 



i=l k=2 ii,...,ife=l 

d n—1 d 

■,ik\ 



< Ci|<^,|V sup \Bl\+C2y] Y] sup IV't'' 

0<t<T , — ' . , 0<t<T 
1=1 — — k=l ii,...,ik=l 

Thus by Gronwall's lemma, we have 

(n-l d \ ( 

V V sup I^Ar'-''"! exp<^ciV sup \Bl 
. ^ , 0<i<T / ^ 0<t<T 

k=l ii,...,«fe=l / k «=1 

This inequality holds true for all < s < L Thus 

(n-l d \ { 

V V sup l^/^r'-'*"! exp<^ciV sup \B\\ 
. , 0<t<T I , 0<t<T 

k=l ti,...,ik=l / k i=l 

which implies (H3l) . □ 



4.2. Malliavin derivative. This subsection is devoted to enhance our Proposition 14. 3[ 
and prove that the Malliavin derivative of yt has also bounded moments of any order in 
our particular nilpotent situation. Notice once again that the boundedness of moments 
of the Malliavin derivative is still an open problem for rough differential equations in the 
general case. We refer to Section [2.3.21 for notations on Malliavin calculus. 

Theorem 4.4. Let the vector fields Vi, 1 = 1,2,--- ,d he smooth with all derivatives 
hounded, satisfying Hypothesis 4jJ_ Assume that all the Lie hrackets of order greater or 
equal to 2 are constants. Then the Malliavin derivative yt has moments of any order. 
More precisely, for any q > 1 and T G (0, oo), 



E 



sup \D^yt\ 

0<u<t<T 



< oo . (45) 



Proof. Go back to our representation f l44j) . which can easily be differentiated in the Malli- 
avin calculus sense in order to obtain 

d n—1 d 

dsD^<f), = 5^vr,(0,)5,*D„0, + J] Yl v^r'-"'=vv,,,...,,,(0,)D„0, 

i=l k=2 ii,...,ife=l 

d n—1 d 

i=l k=2 ji,.--,ifc=l 

where we have set VVj^ , . for the (matrix valued) gradient of V^^ j^, and where we 
recall that the notation has been introduced at Section [2.3.21 
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Since we assume that all the Lie brackets of order greater or equal to 2 of the vector 
fields Vi are constant vector fields, it is easily checked that 

d d n—1 d 

k=2 n,...,ifc=l 



i=l 



Therefore there exist two positive constants Ci,C2 such that 

d 



(46) 



< cij2\Bi\\Du<ps\ 

i=l 

d 



n-l 



Ik \'rsj 



1=1 



k=2 ii,...,jfc=l 



By Gronwall's lemma we obtain 



sup 

o<s<i.o<M<t<r 



\Du(j)s\ < C2exp S sup 15^*1 > X <^ 



i=l 
n—1 d 



0<t<T 



i=l 



1 + sup 

0<s<l 



X 1+ sup V 



0<t<T 



k=2 n,...,jfc=l 



sup 

0<s<l,0<u<t<T 



IV 



i\,...,tk\'Vs) 



Thus 



sup \Duyt\ < C2exp < ci ^ sup |5. 



0<u<t<T 



■1=1 



0<t<T 



X 



E 

i=l 



1+ sup \yt\ 

0<t<T 



n-l 



k=2 ii,...,ik = l 



sup 

0<s<l,0<u<t<T 



tl,...,lk 

t 



□ 



which ends the proof easily by boundedness of moments for yt, Dui^t''"' ^^'^ 



Example 4.5. A classical example of nilpotent vector fields in M"^ is due to Yamato 
Let us check that this example fuUfils our standing assumptions. Indeed, the example 
provided in [31] is the following: 



d d 

= , ^2 = — + 2X2 



dxi 



and A'i 



d ^ d 
- 2xi- 



dx 



dx:> 



Then 



-4- 



dx^ 



[[A,,As],A2] = [[A2,As],A;] = 0. 



It is thus readily checked that the conditions of Theorem 14.41 are met for these vector 
fields. Moreover, in this particular case the solution to equation 022 p is explicit and we 
have 

yl = 1/1 + Bl y^ = y, + Bf, yl = y, + 2 {^If - BoY') , 
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if the solution starts from the initial condition (yi,2/2;2/3)- Interestingly enough, though 
the solution is explicit here, the smoothness of the density of yt is not immediate and we 
recover here the results of j9]. 

4.3. Stochastic flows. The probabilistic proof of the smoothness of density for diffusion 
processes originally given by Malliavin j23] heavily relies on stochastic flows methods and 
their relationship with stochastic derivatives. We now establish those relations for SDEs 
driven by a fractional Brownian motion. 

To this aim, denote by y*''^ the solution to equation (122|) starting from the initial 
condition Us = a aX, time s: 

d 

dyr = Y,yriytldBl te[s,r], yr = a. (47) 

i=l 

The above equation gives rise to a family of smooth nonlinear mappings $^,4 '■ IR*" — ^ W^, 
< s < t <T, determined by $st(a) := y^'"', and the family {^s,u ^ s <t < T} has the 
following flow property (we refer e.g. to [12] for the properties of flows driven by rough 
paths quoted below): 

^s,t = ^u,t O <^s,u , \/0<s<u<t<T. 

Let Js^t denote the gradient of the nonlinear mapping ^s,t with respect to the initial 
condition. Then the family {Js,t', < s < t < T} also satisfies the relation Js^t = Ju,tJs,u 
for < s < M < t < T. In addition, the map Js.t is invertible, and we have Jg t = Jq i-JqI- 
The equation followed by Jq i is obtained by differentiating formally equation (H7|) with 
respect to the initial value a, which yields 

d 

dJo,t = ^Vi{yt)Jo,t dBi , Jo,o = / • 

i=l 

By applying the rules of differential calculus for rough paths, we also get that J^l is 
solution to the following equation: 

d 

dJ,i = -Y.VV.{yt)J^}dBl, J,~i = I. (48) 

i=l 

We have thus ended up with two linear equations for the derivatives of the flow. In our 
nilpotent case, we are able thus able to bound these derivatives along the same lines as 
for Theorem 14.41 

Theorem 4.6. Let the vector fields Vi, 1 = 1,2, he smooth with all derivatives 

hounded, satisfying Hypothesis \4-l\ Assume that all the Lie hrackets of order greater or 
equal to 2 are constant. Then the Jacohian Jq^i md its inverse Jq] have moments of any 
order: for any q>l and T G (0, oo), 



E 



sup |Jo,f 

0<t<T 



< oo, and E 



I 7-1 I "? 

sup Jo,t 
0<t<T 



< OO . (49) 



Proof. As mentioned in the proof of Proposition 14. 3^ one can write $o,f('2) = exp(Z()(a) = 
01 (a), where 0i(a) satisfies (144|) . Thus if we introduce Js = V0s, then Jo^t = Ji where Js 
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satisfies 

d n—1 d 

d,l = J2BlVV,{<Ps)Js + J2 E V'^'-^'^ VV,„...,,,(0,) J,, Jo = /. 

i=l k=2 ji,...,ifc=l 

The first part of (H^ is tlius proved following tlie steps of Proposition 14.31 

As far as tlie second part of f H9|) is concerned, observe tliat JqI = Ji, wliere Jg is 
inverse of Jg. It is clear that satisfies 

d n—1 d 

dsJs = -Y,'^v,{(Ps)JsBi-j2 Yl v^r'-'^'VY,,,,,,,,, J,, jo = /. 

i=l fc=2 ii,...,ifc=l 

Once again, the methodology of Proposition 14.31 easily yields our claim. 

□ 

Corollary 4.7. Under the same assumptions as in Theorem \4.6[ the following holds true: 

(i) For any < 7 < and q > 1 we have 

E [\\Dy,r^J + E Ol^o.ll^] + E [ll^^.l^ J < ct,„ (50) 
for a finite constant CT,q- 

(ii) As a consequence, inequality ([J^j also holds true when the \\ ■ ||^^oo norms are replaced 
by norms in Ti, where T-i has been defined at Section \2.3.1[ 

(iii) For any smooth bounded vector field U onW^ andt E [0,T], set Z^^ := {jQlU{yt), rj). 
Then is a controlled process, and satisfies the inequality E[Af'^[Z^] Q'^{W^)]] < CT,q 
for any q > 1 and a finite constant CT,q- 

Proof. Going back to equation ( 146|) . it is readily checked that all the terms u ^ Duipl^''"'^'' 
are C^-Holder continuous on [0,t] for any < H, since the elements ^jj^'---'^^ are nice 
multiple integrals with respect to B. Moreover, we have 

E J < 00, 

for any m > 1. This easily yields E[||Z}y^^''"'*'°||^ < 00 by a standard application of 
Gronwall's lemma, as in the proof of Theorem 14.41 

Our second assertion stems from the fact that one can choose 1/2 — H < < H, since 
H > 1/3. For such 7 we have Cj C H, which ends the proof. 

Finally, our claim (iii) derives from the fact that the equation governing is of the 
following form: 




The process Z can thus be decomposed as a controlled process as in Section \2.2\ and 
since we already have estimates for JqI and yt, the bound on 'E[J\f'^[Z'^; Q'^(]R™)]] follows 
easily. 

□ 
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4.4. Proof of Theorem 11.11 As mentioned in the introduction, once we have shown 
the integrabihty of the MaUiavin derivative and proved a Norris type lemma, the proof of 
our main theorem goes along classical lines. We have chosen to follow here the exposition 
of [16], to which we refer for further details. 

Step 1: Reduction to a lower hound on Holder norms. Recall that the process has 
been defined for any smooth vector field U in Corollarv 14.71 For any p > 1, our first goal 
is to reduce our problem to the existence of a constant Cp such that 

P(||^''NU,oo<£) <Ope^\ (52) 

for 1 < A; < (i, a given a G (1/3, 1/2), all e G (0, 1) and where we observe that all the 
norms below are understood as norms on [0,T]. 

Indeed, according to [T6| Relation (4.9)] the smoothness of density can be obtained 
from the estimate 

where we recall that "H has been defined at Section 12.3.11 Furthermore, we have 

P <e)<V [WZ^'-W,. <e)<V (llZ^'^IUi < e) . 

It is thus sufficient for our purposes to check 

P(||Z^^|Ui <£) <Cp£^'. (53) 

In order to go from (153|) to ( 152|) . let us use our interpolation bound (123|) in the following 
form: for any < ?7 < 1 and 0<a;<p<-frwe have 

II&IIli > V/^^~"Ull&IU,oo-C„,,r/||6||,,,o). 
Take now 5 G (0, 1) to be fixed later on and r/^(^'~") = e^"^, that is 77 = Then 

P {WZ^Hl^ < e) < P (||^''1U,oo < 2£^) + R, where /2 = P (^||Z^M|p,oc < , 

(54) 

with V = p — a — {l — {p — a))6. Choose now 5 small enough, so that u > 0. Since 
admits moments of any order according to Corollary 14. 7^ it is easily checked that R can 
be made smaller than any quantity of the form CgS'^. It is thus sufficient to prove (1521) in 
order to get the smoothness of density for yt. 

Step 2: An iterative procedure. For / > 1 and x G M™, let V;(x) be the vector space 
generated by the Lie brackets of order / of our vector fields Vi, . . . , at point x: 

Vi{x) = Span {[\4, ■ ■ ■ \4j(x); ] < lA < . . . < d] . 

We assume that the vector fields are £-hypoelliptic for a given £ > 0, which can be read 
as Vi{x) = for any x G M™. In order to start our induction procedure, we set ai = a, 
so that we have to prove PdlZ^'' ||q,^^oo < £^) < CpE^. 
Recall that Z^'= satisfies the relation 

Zr^ = {r^,V,{a)) + f2 f ZY^-^^^dBl 

Thus Proposition 13.21 asserts that for any 1/3 < a2 < ai < H there exists q2 > such 
that 
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Relation fl52p is thus implied by 

Iterating this procedure we end up with the following claim to prove: Bii[e) < CpE^ for all 
e e (0,1), with 

where the intersection above extends to all possible combinations 1 < ki, . . . , < d, and 
where 1/3 < ai < ■ ■ ■ < ai < H. 

Going back now to the very definition of as Z^^ = {jQlU{yt), 77), it is readily checked 
that 

B,{e) < P {{v, VM) <e^...Av. [Vk, ■ ■ ■ V^) < e^O • 

Owing to the fact that V£{a) = M'", we thus have Bi{e) = for e small enough, which 
ends the proof. 

Acknowledgement: We would like to thank Ivan Nourdin for pointing us out the use 
of Hermite polynomials in the computation of 4*^^ order variations (see proof of Proposi- 
tion [32] Step 2). 
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